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A spatially continuous force-based model for simulating pedestrian dynamics is introduced which 
includes an elliptical volume exclusion of pedestrians. We discuss the phenomena of oscillations 
and overlapping which occur for certain choices of the forces. The main intention of this work is 
the quantitative description of pedestrian movement in several geometries. Measurements of the 
fundamental diagram in narrow and wide corridors are performed. The results of the proposed 
model show good agreement with empirical data obtained in controlled experiments. 



I. INTRODUCTION 

For a beneficial application of pedestrians dynamics, 
robust and quantitatively verified models are required. 
A wide spectrum of models has been designed to sim- 
ulate pedestrian dynamics. Generally these models can 
be classified into macroscopic and microscopic models. 
In macroscopic models the system is described by mean 
values of characteristics of pedestrian streams e.g., den- 
sity and velocity, whereas microscopic models consider 
the movement of individual persons separately. Micro- 
scopic models can be subdivided into several classes e.g., 
rule-based and force-based models. For a detailed discus- 
sion, we refer to [TJ [5] . In this work we focus on spatially 
continuous force-based models. 

Force-based models take Newton's second law of dy- 
namics as a guiding principle. Given a pedestrian i with 

coordinates Ri we define the set of all pedestrians that 
influence pedestrian i at a certain moment as 

M := {j : || - % ||< r c A i "feels" j} (1) 

where r c is a cutoff radius. We say pedestrian i "feels" 
pedestrian j if the line joining their centers of mass does 
not intersect any obstacle. In a similar way we define the 
set of walls or borders that act on pedestrian i as 

W l :={w:\\lC t -Tt l \\<r c } (2) 

where lOj £ w is the nearest point on the wall w to the 
pedestrian i. 

Thus, the movement of each pedestrian is defined by 
the equation of motion 

m i t = f i = Ff^+Y l ^?+ E ( 3 ) 



where denotes the repulsive force from pedestrian j 
acting on pedestrian i, is the repulsive force emerg- 
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ing from the obstacle w and Ff rv is a driving force, m, 
is the mass of pedestrian i. 

The repulsive forces model the collision-avoidance per- 
formed by pedestrians and should guarantee a certain 
volume exclusion for each pedestrian. The driving force, 
on the other hand, models the intention of a pedestrian 
to move to some destination and walk with a certain de- 
sired speed. The set of equations |3]) for all pedestrians 
results in a high-dimensional system of second order or- 
dinary differential equations. The time evolution of the 
positions and velocities of all pedestrians is obtained by 
numerical integration. 

Most force-based models describe the movement of 
pedestrians qualitatively well. Collective phenomena like 
lane formation [3H5] , oscillations at bottlenecks [3J 0] , the 
"faster-is-slower" effect [5J [7] , clogging at exit doors [H [5] 
are reproduced. These achievements indicate that these 
models are promising candidates for realistic simulations. 
However, a qualitative description is not sufficient if reli- 
able statements about critical processes, e.g., emergency 
egress, are required. Moreover, implementations of mod- 
els often require additional elements to guarantee realis- 
tic behavior, especially in high density situations. Here 
strong overlapping of pedestrians [5J |5] or negative and 
high velocities [3j [8] occur which then has to be rectified 
by replacing the equation of motion ([3| by other proce- 
dures. 

Force-based models contain free parameters that can 
be adequately calibrated to achieve a good quantitative 
description [9Hl3|. However, depending on the simulated 
geometry the set of parameters often changes. In most 
works quantitative investigations of pedestrian dynamics 
were restricted to a specific scenario or geometry, like 
one-dimensional motion [T3J j behavior at bottlenecks [TTJ 
HHJ HO]; two-dimensional motion 12J or outflow from a 
room [TfjHTg] . 

In this work we restrict ourselves to corridors and ad- 
dress the possibility of describing the movement of pedes- 
trians in wide and narrow corridors reasonably and in a 
quantitative manner with a unique set of parameters. At 
the same time, the modelling approach should be as sim- 
ple as possible. 
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FIG. 1. (Color online) Direction of the repulsive force. 



In the next section, we propose such a model which 
is solely based on the equation of motion (j3J). Further- 
more the model incorporates free parameters which allow 
calibration to fit quantitative data. 



affected by the presence of faster pedestrians in front of 
them: 



J ( Vi - Vj ) ■ etj if (vi - Vj) ■ eij > 







otherwise. 



(7) 



As in general pedestrians react only to obstacles and 
pedestrians that are within their perception, the reac- 
tion field of the repulsive force is reduced to the angle 
of vision (180°) of each pedestrian, by introducing the 
coefficient 



eij)/vi 



if v, 



otherwise. 



-ij 



> & Vi . + 



(8) 



II. THE CENTRIFUGAL FORCE MODEL 

The Centrifugal Force Model (CFM) [5] takes into ac- 
count the distance between pedestrians as well as their 
relative velocities. Pedestrians are modelled as circular 
disks with constant radius. Their movement is a direct 
result of superposition of repulsive and driving forces act- 
ing on the center of each pedestrian. Repulsive forces 
acting on pedestrian i from other pedestrians in their 
neighborhood and eventually from e.g. walls and stairs 
to prevent collisions and overlapping. The driving force, 
however, adds a positive term to the resulting force, to 
enable movement of pedestrian i in a certain direction 
with a given desired speed v® . The mathematical expres- 
sion for the driving force is given by 



-idrv 



(4) 



with a time constant r. 

Given the direction connecting the positions of pedes- 
trians i and j: 



Rij Rj Ri ) 



!?<, 
Ri 



The repulsive force then reads (see Fig. [IJ 



n t j 



(5) 



(6) 



This definition of the repulsive force in the CFM reflects 
several aspects. First, the force between two pedestrians 
decreases with increasing distance. In the CFM it is in- 
versely proportional to their distance Rij. Furthermore, 
the repulsive force takes into account the relative velocity 
Vij between pedestrian i and pedestrian j. The following 
special definition provides that slower pedestrians are not 



The coefficient kij is maximal when pedestrian j is in the 
direction of movement of pedestrian i and minimal when 
the angle between j and i is bigger than 90°. Thus the 
strength of the repulsive force depends on the angle. 

As mentioned earlier the CFM is complemented with 
a "Collision Detection Technique" (CDT) to manage 
conflicts and mitigate overlappings between pedestrians. 
Fig. [2] depicts schematically the definition of the CDT. 
Although CDT is relatively simple, it adds an amount of 
complexity to the initial model defined with Eq. ^ and 
masks the main idea behind the repulsive forces. In the 
following we systematically modify the expression of the 
repulsive force to enable a better quantitative description 
of pedestrian dynamics. 



III. OVERLAPPING VS. OSCILLATION 

In this work we consider a velocity-dependent vol- 
ume exclusion of pedestrians. Overlapping between two 
pedestrians occurs when their geometrical form (circle, 
ellipse, ...) overlaps. Modelling a pedestrian as a circle 
or ellipse is just an approximation of the human body. 
Therefore, a certain amount of overlapping could be ac- 
ceptable and might be interpreted as "elastic deforma- 
tion". However, for the deformed body the center of 
mass no longer coincides with the center of the circle or 
ellipse. For this reason overlapping is a serious problem 
that should be dealt with. 

In [5T] it was shown that the introduction of a CDT 
is necessary to mitigate overlapping among pedestrians. 
The CDT keeps pedestrians away from each other with 
a distance of at least r, where r represents the radius of 
the circle modelling the volume exclusion of pedestrians. 

Our goal is to simplify the model by dispensing with 
the CDT and improve the repulsive force to compensate 
for the effects of the missing CDT on the dynamics. To 
introduce the shape of the modeled pedestrians in Eq. ^ 
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FIG. 2. Schematic representation of the collision detec- 
tion technique (CDT), which is an important component in 
the CFM [S], to manage collisions and mitigate overlapping 
among pedestrians. In our model we do not need the CDT, 
which is a considerable simplification in comparison to the 

CFM pg. 



we transform the singularity of the repulsive force from 
to 2r: 



= —rriik. 



Rij ~ 2r 



(9) 



Due to the quotient in Eq. ([9| when the distance is small, 
low relative velocities lead to an unacceptably small force. 
Consequently, partial or total overlapping is not pre- 
vented. Introducing the intended speed in the numera- 
tor of the repulsive force eliminates this side-effect. This 
dependence on the desired speed is motivated by the ob- 
servation that for faster pedestrians stronger repulsive 
forces are required to avoid collisions with other pedes- 
trians and obstacles. Thus, the repulsive force is changed 
to 



F p = 



R, 



2r 



(10) 



with a free parameter rj to adjust the strength of the 
force. 

Those two changes in the repulsive force cause the 
emergence of two phenomena: Overlapping and oscilla- 
tions. In the following we will define quantities to study 
those phenomena. 

Avoiding overlapping between pedestrians and oscil- 
lations in their trajectories is difficult to accomplish in 




FIG. 3. The overlapping area between pedestrians i and j 
varies between and 1. 



force-based models. On one hand, increasing the strength 
of the repulsive force with the aim of excluding overlap- 
ping during simulations leads to oscillations in the trajec- 
tories of pedestrians. Consequently backward movements 
occur, which is not realistic especially in evacuation sce- 
narios. 

On the other hand, reducing the strength of the re- 
pulsive force (to avoid oscillations) leads inevitably to 
overlapping between pedestrians or between pedestrians 
and obstacles. 

To solve this dilemma one has to find an adequate value 
of the strength of the repulsive force: it should neither 
be too high so that oscillations will appear, nor too low 
so that overlapping will be observed. 

To understand this duality we quantify overlapping 
and oscillations during simulations. First, we define an 
overlapping-proportion during a simulation as: 



» _ 



^ t=t ond i=N j=N 

- E E E °« > 

t—0 2—1 j>t 



with 



A; 



mm(j4j, Aj) 



< 1, 



(11) 



(12) 



where N is the number of simulated pedestrians. Aij is 
the overlapping area of the circles i and j with areas Ai 
and Aj , respectively (see Fig. |3| . n ov is the cardinality 
of the set 



O : K, : <>,, £ 0} . 



(13) 



For n ov = 0, o^ v ' is set to zero. 

For a pedestrian with velocity ~Vi and desired velocity 

we define the oscillation-proportion as 

t=t cnd j=JV 



L OS . n -i 
t—0 2—1 



(14) 



where Si quantifies the oscillation-strength of pedestrian 
i and is defined as follows: 



(15) 
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FIG. 4. (Color online) Oscillation-proportion o ( - 3 ^ and 
overlapping-proportion o w as function of the interaction 
strength r\ obtained from 200 simulations with different initial 
conditions. Oscillations increase with increasing strength of 
the repulsive force, while overlaps become negligible for larger 
77. The case rj — is the CFM. In each run the simulations 
for different rj are started with the same initial values. 



with 



Vi ■ Vi 

and n os is the cardinality of the set 
S := \s, : fli £ 0}. 



(16) 



(17) 



Here again is set to zero if n os = 0. The propor- 
tions o^) and are normalized to 1 and describe the 
evolution of the phenomena overlapping and oscillations 
during a simulation. 

In order to exemplify the behavior of these two coupled 
phenomena we simulate an evacuation of 35 pedestrian 
from a 4m x 4m room with an exit of 1.2m and deter- 
mine o^ v ' and for different values of rj in Eq. (10 1 



Results are shown in Fig. |3J 77 = is a special case of 
the model and represents the CFM [5]. The high values 
of the overlapping proportion suggest that simulations 
using only CFM without the CDT lead to unreasonable 
results. For further details we refer to |21) . 

The introduction of the intended velocity in the repul- 
sive force enhances the ability of the repulsive force to 
guarantee the volume exclusion of pedestrians. This is 
reflected by the decreasing of th e overlapping-proportion 
o^"> while increasing 77 (Eq.[l0j|. See Fig. [3] 

Meanwhile, the oscillation-proportion increases, 
thus the system tends to become instable. Large val- 
ues of the oscillation-proportion o^ 3 ' imply less stability. 
For Si — 1 one has = —vt°, i.e. a pedestrian moves 
backwards with desired velocity. Even values of Si higher 



than 1 are not excluded and can occur during a simula- 
tion. Therefore, a careful calibration of 77 is required to 
achieve an optimal balance between overlapping and os- 
cillations. 

Unfortunately, it is not possible to adjust the strength 
of the repulsive force by means of rj in order to get an 
overlapping-free and meanwhile an oscillation-free simu- 
lation. Nevertheless, by proper choice of 77 one can reduce 
the amount of overlapping among pedestrians such that 
it becomes negligible and can be interpreted as a defor- 
mation. This characteristic of the GCFM is not fulfilled 
by the CFM [5], where total overlapping (o^ = 1) can be 
observed. 

Furthermore, the quantities and o^"' provide a cri- 
terion to choose an optimal value for 77, which is given by 
the intersection of the curves representing and o^ v \ 



IV. HARD CIRCLES VS. DYNAMICAL 
CIRCLES: THE FUNDAMENTAL DIAGRAM 
FOR SINGLE FILE MOVEMENT 

It is suggested that the effective space requirement 
of a moving pedestrian varies with velocity. Usually, 
the projection of the pedestrian's shape to the two- 
dimensional plane is modeled as a circle with a radius 
r [3l [3 [10] . Thompson suggested a three-circle represen- 
tation for main body and shoulders [22:. According to 
|23j . however, the radius of the circle varies such that the 
space requirement of pedestrians increases significantly 
as speed increases. In [14) a linear velocity-dependence 



T r Vi 



(18) 



of the radius with parameters r m i n and r r was sug- 
gested. "Space requirement" encompasses the physical 
area taken by the torso together with the motion of the 
legs, lateral swaying and a safety margin. 
The repulsive force reads 



F; 



with 



dist; 



TTli hi 



dist, 



n(vi) -r 3 (vj) 



(19) 



(20) 



the effective distance between pedestrian i and j and n 
the radius of pedestrian i as defined in Eq. (18 1. 



V. ELLIPTICAL VOLUME EXCLUSION OF 
PEDESTRIANS 

One drawback of circles that impact negatively the 
dynamics is their rotational symmetry with respect to 
their centers. Therefore, they occupy the same amount 
of space in all directions. In single file movement this is 
irrelevant since the circles are projected to lines and only 



the required space in movement direction matters. How- 
ever, in two dimensional movement the aforementioned 
symmetry lasts by occupying unnecessary lateral space. 

In [21] Fruin introduced the "body ellipse" to describe 
the plane view of the average adult male human body. 
Pauls presented ideas about an extension of Fruin's 
ellipse model to better understand and model pedestrian 
movement as density increases. Templer |25| noticed that 
the so called "sensory zone" , which is a bubble of space 
between pedestrians and other objects in the environ- 
ment to avoid physical conflicts and for psychocultural 
reasons, varies in size and takes the shape of an ellipse. 
In fact, ellipses are closer to the projection of required 
space of the human body on the plane, including the ex- 
tent of the legs during motion and the lateral swaying of 
the body. 

Having the ambition to describe with the same set of 
parameters the dynamics in one- and two-dimensional 
space we extend our model by introducing an elliptical 
volume exclusion of pedestrians. Given a pedestrian i we 
define an ellipse with center (xi,yi), major semi-axis a 
and minor semi-axis b. a models the space requirement 



in the direction of movement. In analogy to Eq. (18) we 

set 



+ T a Vi 



with two parameters a m i n and T a . 

Fruin [24] observed body swaying during both human 
locomotion and while standing. Pauls |26| remarks that 
swaying laterally should be considered while determining 
the required width of exit stairways. In [20] character- 
istics of lateral swaying are determined experimentally. 
Observations of experimental trajectories in |20j indicate 
that the amplitude of lateral swaying varies from a max- 
imum 6 max for slow movement and gradually decreases 
to a minimum 6 m i n for free movement when pedestrians 
move with their free velocity (Fig. [5| . Thus we describe 
with b the lateral swaying of pedestrians and set 



b = K 



(^max frmin) q 



(22) 



Since a and b are velocity-dependent, the inequality 



b < a 



(23) 



does not always hold for the ellipse i. In the rest of this 
work we denote the semi-axis in the movement direction 
by a and its orthogonal semi-axis by b. 



VI. ELLIPTICAL VOLUME EXCLUSION AND 
FORCE IMPLEMENTATION 

In this section we give some mathematical insights con- 
cerning the implementation of the repulsive forces. 




FIG. 5. Off-line trajectory detection with PeTrack [27] . Left: 
The trajectory of the detected pedestrian shows strong sway- 
ing. Right: The faster pedestrians move, the smoother and 
weaker is the swaying of their trajectories. 



A. Repulsive Forces between Pedestrians 

In order to calculate the repulsive force emerging from 
pedestrian j acting on pedestrian i according to Eq. ( 19 ) 
we require the distance between the borders of the el- 



lipses, along a line connecting the two pedestrians dist^. 
See App.f^jfor more details on disty. 

Another important quantity is the distance of closest 
approach or contact distance of two ellipses / which is the 
(21) minimum of disty while i and j are not overlapping. Un- 
like for circles, I can be non-zero for ellipses and depends 
on their orientations. In [28j an analytical expression for 
the distance of the closest approach of two ellipses with 
arbitrary orientation is derived. Fig. [7] shows how distij 
and I goes in the repulsive force. 



B. Repulsive Forces between Pedestrians and Walls 

The repulsive force between a pedestrian i and a wall 
is zero if i performs a parallel motion to the wall. While 
this behavior of the force is correct, it leads to very small 
repulsive forces when the pedestrians motion is almost 
parallel to the wall. For this reason we characterize in this 
model walls by three point masses acting on pedestrians 
within a certain interaction range (Fig. |6| . The middle 
point is the point with the shortest distance from the 
center of the pedestrian to the line segment of the wall. 
All three points have to be computed at each step as 
the pedestrian moves. The distance between the three 
wall points is set to the minor semi-axis of an ellipse If 
one lateral point (wi + i or lOj+i ) does not lie on the line 
segment of the wall, then it will not be considered in the 
computation of the repulsive force. 

The number of point masses have been chosen by a 
process of trial and error. Simulations have shown that 
three point masses are sufficient to keep pedestrians away 
from walls. Meanwhile they are computationally cost- 
effective. 

As walls are static objects, the repulsive force emerging 



G 




s / r eps r c r c 



from a wall w and acting on pedestrian i simplifies to 



i+l 

£ 

.3=1-1 



(24) 



with 



Fi 



dist,; 



-e^, j E {i-1, i, i+l} . 

(25) 

v " is the component of the velocity normal to the wall, 
ki Wi and as defined resp. in Eqs. (Tsl) and |5~1) in 
Sec. HE 

The distance between a line w and the ellipse i is 



dist. 



(26) 



with Ti the polar radius determined in Eq. ( A3 1 and k 



the distance of point Oi to the line w. Further details can 
be found in App. [b] According to the distance I defined 
for the repulsive forces between pedestrian in Sec. |VI A| 
we introduce the distance of the closest approach between 
an ellipse and a line k, see App. [b] and Fig. [7] for details. 
Note that in Eq. ( 25 ) , ki Wi in the force is independent 



of the chosen lateral wall point Wj. That means, if a 
pedestrian is moving parallel to the wall, ki W . — and 
thus the points j — 1 and j + 1 have no effects. 



C. Numerical Stabilization of the Repulsive Force 

In this section we describe a numerical treatment of 
the repulsive force. For the sake of simplicity, we focus 
on the case of pedestrian-pedestrian interactions. The 
pedestrian-wall case is treated similarly. 

The strength of the repulsive force decreases with in- 
creasing distance between two pedestrians. Nevertheless 
the range of the repulsive force is infinite. This is unreal- 
istic for interactions between pedestrians. Therefore, we 
introduce a cut-off radius r c — 2 m for the force limiting 
the interactions to adjacent pedestrians solely. To guar- 
antee robust numerical integration a two-sided Hermitc- 
interpolation of the repulsive force is implemented. The 
interpolation guarantees that the norm of the repulsive 



FIG. 7. (Color online) The interpolation of the repulsive force 



between pedestrians i and j Eq. ( 19 1 dep ending on distij and 



the distance of closest approach I, see Sec. 



VIA 



As the repul- 
sive force also depends on the relative velocity Vij, this figure 
depicts the curve of the force for wy = const. The left and 
right dashed curves are defined in Eqs. ( |28[ ) and (27} respec- 
tively. The wall-pedestrian interaction has an analogous form 
with distij and I replaced by distjj and k, respectively. 



force decreases smoothly to zero for disti 



For 



distjj — > l + the interpolation avoids an increase of the 
force to infinity but to f m — 3F , 4 rcp (r cps ) at s 



' cps 



and r. 



cps 



0.1 m, where it remains constant, distjj and I 
are illustrated in Sec. | VI A} Fig. [7] shows the dependence 
of the repulsive force on the distance for constant relative 
velocity. 

The right interpolation function P r and the left one Pi 
(dashed parts of the function in Fig. [7]) are defined using 



P r (r c )= F^ p (f c ), P r {r c ) = 
(PrY(fc)= (F^)' (f c ), (P r )'(r c ) = (27) 



with r c = r c — r cps and 

Pi{so)= fm, Pi(r cps ) = F-- P (r cps ) 
(P*)'(4)= 1, (Pi)'(r eps ) = (F^Y (r cps ) . (28) 

where the prime indicates the derivative, sq is the min- 
imum allowed magnitude of the effective distance of two 
ellipses. Due to the superposition of the forces the in- 
equality: 

dist^ > s . (29) 
for pedestrians i and j is not guaranteed. 



VII. SIMULATION RESULTS 

The initial value problem in Eq. ([3]) was solved using an 
Euler scheme with fixed-step size At = 0.01 s. First the 
state variables of all pedestrians are determined. Then 
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FIG. 8. (Color online) Top: Velocity-density relation for one- 
dimensional movement compared to experimental data [29] . 
For the simulations, r a is set to 0.53 s . Bottom: Changing 



T a in Eq. (21 1 influences the slope of the diagram. a m i n has 
been kept equal to 0.18 m. r a = represents pedestrians with 
constant space-requirement. 



the update to the next step is performed. Thus, the 
update in each step is parallel. 

The desired speeds of pedestrians are Gaussian dis- 
tributed with mean fi = 1.34 m/s and standard deviation 
a = 0.26 m/s. The time constant r in the driving force 
Eq. Q is set to 0.5 s, i.e. r ^> At. For simplicity, the 
mass rrii is set to unity. 

In order to verify the model and evaluate the differ- 
ence of the elliptical shape of the volume exclusion ver- 
sus the circular one we measure the fundamental diagram 
in two-dimensional space with the same set of parame- 
ter as for the one-dimensional fundamental diagram. In 
the one-dimensional case only the space requirement of 
pedestrians in movement direction, expressed in terms of 
the semi-axis a, influences the dynamics of the system. 
We set a m i n 



FIG. 9. (Color online) Density-velocity relation in a corridor 
of dimensions 25 m x 1 m in comparison with experimental 
data obtained in the HERMES-project [30]. For the simula- 
tion with circles, b is set to be equal to a. 



To illustrate the impact of the velocity-dependence of 
the radius on the dynamics of pedestrians we measure 
the one-dimensional fundamental diagram in a corridor 
of 26 m with periodic boundary conditions. The mea- 
surement segment is 2 m long and situated in the middle 
of the corridor. Details about the measurement method 
are given in App. [C] 

The results for the one-dimensional fundamental dia- 
gram are shown in Fig. [8] and compare well with experi- 
mental data. Ellipses with velocity-dependent semi-axes 
emulate the space requirement of the projected shape of 
pedestrians better. Even the shape of the fundamen- 
tal diagram is reproduced after inclusion of this velocity- 
dependence. 

We extend the simulation to two-dimensional space 
and simulate a 25 m x 1 m corridor with periodic bound- 
ary conditions. A measurement segment of 2 m x 1 m was 
set in the middle of the corridor. We use the same mea- 
surement method as for the single-file case (see App. [C|. 
Calibration of the parameters of the lateral semi-axis b 
(&min and 6 ma x in Eq. 22 ) leads to the values 6 m ; n = 0.2 m 
and b max = 0.25 m. The simulation result is shown in 
Fig [9} 



0.18 m and r a = 0.53 s (see Eq. 21) 



With the chosen dimensions of the semi-axes a and b 
the model yields the right relation between velocity and 
density both in single-file movement and wide corridors, 
although only a corridor width of 1 m was investigated. 
One remarks that the fundamental diagram for elliptical 
shaped particles is an upper bound for that of circular 
ones, especially at low and medium densities. At high 
densities there is no noticeable difference between both 
shapes. 



VIII. CONCLUSIONS 

We have proposed modifications of a spatially contin- 
uous force-based model [5] to describe quantitatively the 
movement of pedestrians in one- and two-dimensional 
space. Besides being a remedy for numerical instabili- 
ties in CFM the modifications simplify the approach of 
Yu et al. [5] since we can dispense with their extra "col- 
lision detection technique" without deteriorating perfor- 
mance. The implementation of the model is straightfor- 
ward and does not use any restrictions on the velocity. 
Furthermore, we introduced an elliptical volume exclu- 
sion of pedestrians and studied its influence compared 
to the standard circular one. Simulation results show 
good agreement with experimental data. Nevertheless, 
the model contains free parameters that have to be tuned 
adequately to adapt the model to a given scenario. Fur- 
ther improvement of the model could be made by includ- 
ing, for example, a density-dependent repulsive force. 

Although the model describes quantitatively well the 
operative level of human behavior, it does not consider 
aspects of the tactical and strategic levels [31] , Phenom- 
ena like cooperation, changing lanes and overtaking are 
not reproduced, especially in bi-directional flow. 

The source code of this model will be released under 
the GNU General Public Licence (GPL (32]) and will be 
available for download from |33| . 




FIG. 10. (Color online) disty is the distance between the 
borders of the ellipses i and j along a line connecting their 
centers. 



for the polar radius with 



sin 2 a, 



6? 



(A4) 



In the same manner, we determine the polar radius rj . 
Finally, the distance dist^ between the centers of the 
ellipses i and j is determined as follows (Fig. 10 1: 



disti 



HI O^Oj II 



(A5) 



Note that the distance between two ellipses can be non- 
zero even when the ellipses touch or overlap. 
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Appendix A: Distance between two ellipses 

In this appendix we give details about the calculation 
of the distance dist^ between two ellipses which is defined 
as the distance between the borders of the ellipses, along 
a line connecting their centers (Fig. |l0|). 

By proper choice of the coordinate system the ellipse 
i may be written as quadratic form, 



= 1 . 



(Al) 



In polar coordinates, with the origin at the center of the 
ellipse and with the angular coordinate on measured from 
the major axis, one gets 



(A2) 



x = ri cos(ai) , y — Ti sin(aj) 



By replacing the expressions of x and y in Eq. ( |Al[ ) and 
rearranging we obtain the expression 



1 = 0, 



(A3) 



Distance of closest approach of two ellipses is the small- 
est distance between their borders, along a line connect- 
ing their centers while they are not overlapping. See 
Fig. (11) top. To mitigate overlapping the repulsive 



forces are high for distances in a certain neighborhood 
of the distance of closest approach, see I in Fig. [Tj An 
analytical solution of this distance for two arbitrary el- 
lipses is presented in [2"5] . 

In this appendix we describe an algorithm to calculate 
the distance of closest approach of an ellipse and a line 
(A), which is the distance between the border of the el- 
lipse, along a line connecting its center o and the closest 
point on the line to o. For this purpose consider without 
loss of generality an ellipse i in canonical position and let 
(A') be the line tangential to the ellipse i and parallel to 
(A) (Fig. JTT1) bottom: 



(A): V 



(A'):y = 



d'. 



(Bl) 



with known coefficients c and d. 

To determine d! we solve the intersection equations of 
an ellipse and a line, which yields the quadratic equation 



q'x 2 + p'x + s' 



0, 



with 



c 



led! 



and s 



dT 
b 2 



(B2) 



1. (B3) 
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As (A') is tangential to the ellipse we have 

D = 



(B4) 



with D the discriminant of Eq. (B2|. Solving |B4j) gives 
d' = ±^b 2 + a 2 c 2 . (B5) 





Pit) 



N in (t) 

I m. 



(C3) 



with l m — 2m the length of the measurement area in the 
movement direction and iVjn (i) is the number of pedestri- 
ans within the area at time t. In one dimensional space 
the measurement area is reduced to a measurement seg- 
ment of length l m . 




FIG. 11. (Color online) Top: Distance of closest approach of 
two ellipses. Bottom: Distance of closest approach between 
an ellipse and a line. 



Finally the distance of closest approach of the ellipse i 
and line (A) is 



k — Ti , 

with k\ the distance o f Cj to (A') and r, 
as determined in Eq. (A3). 



(B6) 

the polar radius 



Appendix C: Measurement method 



The mean velocity of pedestrian i that enters the mea- 
surement are at (x^j-yf 1 ) and leaves it at (x° ut ,y° ut ) is 
determined as 



V(x° ut - xf) 2 + (y° ut - yff 

f out _ t in 



(CI) 



where if is the entrance time and t° ut exit time of i. For 
the one-dimensional case yf 1 = y° nt = 0. 
The density is defined as follows: 

Pi = ^-^j^P(t)dt (C2) 
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